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Abstract. 

Understanding axion cosmology has important experimental 
consequences since it constrains the range of allowed values 
for the axion mass. In the standard thermal scenario, which 
assumes Peccei-Quinn symmetry restoration after inflation, an 
axion string network forms at the phase transition T ^ /a and 
then radiatively decays into a cosmological background of ax- 
ions. Under standard assumptions for the evolution of this string 
network and the radiation from it, axions must have a mass 
ma ~ 100 iieV with specified large uncertainties. We discuss 
critically the various suggestions in the literature that the axion 
mass might be lighter. 



1. Introduction 

The axion has remained a popular dark matter candidate because of its 
enduring motivation as an elegant solution to the strong CP-problem Q]. 
Despite early hopes of discovery, it turned out that in order to be consistent 
with accelerator searches and astrophysics, the axion must be nearly 'invis- 
ible' and extremely light. Its couplings and mass are inversely proportional 
to the (large) Peccei-Quinn (PQ) scale /a as in 

„...6^2xl0-ev(i2^). (1) 

Accelerator constraints have been largely superseded by those from astro- 
physics; because the axion is so weakly coupled, volume effects can com- 
pete with other surface and convective stellar energy loss mechanisms. The 
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strongest astrophysical constraints on the axion mass derive from studies 
of supernova 1987a and conservative estimates yield ma ^ 10 meV [Q. The 
present programme of large-scale axion search experiments |^ are sensitive 
to a mass range ttIq ~ 1-10 /xeV, which has been chosen for a variety of 
historical and technological reasons. Here, we discuss the motivation from 
cosmology for a larger mass axion (ma ~ 100/xeV) if it forms a substan- 
tial part of the dark matter of the universe. We will also comment on 
suggestions in the literature that question this point of view. 



2. Standard axion cosmology 

The cosmology of the axion is determined by the two energy scales /a and 
■'^QCD- The first important event is the PQ phase transition which is bro- 
ken at a high temperature T ~ /a<;l0^GeV. This creates the axion, at 
this stage an effectively massless pseudo-Goldstone boson, as well as a net- 
work of axion strings |^ which decays gradually into a background cosmic 
axions d. (Note that one can engineer models in which an inflationary 
epoch interferes with the effects of the PQ phase transition, this possibility 
is discussed in ref. Q). At a much lower temperature T ^ Aqcd after ax- 
ion and string formation, instanton effects 'switch on', the axions acquire a 
small mass, domain walls form Q between the strings Q and the complex 
hybrid network annihilates into axions in about one Hubble time 

There are three possible mechanisms by which axions are produced in 
the 'standard thermal scenario', which assumes PQ symmetry restoration 
after any inflationary epoch: (i) thermal production, (ii) axion string radi- 
ation and (iii) hybrid defect annihilation when T — Age d ■ Axions consis- 
tent with the astrophysical bounds must decouple from thermal equilibrium 
very early; their subsequent history and number density is analogous to the 
decoupled neutrino, except that unlike a lOOeV massive neutrino, thermal 
axions cannot hope to dominate the universe with ma ^ lOmeV. We now 
turn to the two dominant axion production mechanisms, but first we ad- 
dress an important historical digression. 

2.1. Misalignment misconceptions 

The original papers on axions suggested that axion production primarily 
occurred, not through the above mechanisms, but instead by 'misalignment' 
effects at the QCD phase transition Q. Before the axion mass 'switches 
on', the axion field 6 takes random values throughout space in the range 
to 27r; it is the phase of the PQ-field lying at the bottom of a U{1) 'Mexican 
hat' potential. However, afterwards the potential becomes tilted and the 
true minimum becomes = 0, so the field in the 'misalignment' picture 
begins to coherently oscillate about this minimum; this homogeneous mode 
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corresponds to the 'creation' of zero momentum axions. Given an initial 
rms value 9i for these oscillations, it is relatively straightforward to estimate 
the total energy density in zero momentum axions and compare these to 
the present mass density of the universe (assuming a flat = 1 FRW 



where A w 3 accounts for both model-dependent axion uncertainties 
and those due to the nature of the QCD phase transition, and h is the 
rescaled Hubble parameter. The function f{9) is an anharmonic correction 
for fields near the top of the potential close to unstable equilibrium « tt, 
that is, with f(0) = 1 at the base 9 and diverging logarithmically for 
^ TT 1^. If valid, the estimate (^) would imply a constraint TOa^5/ieV 
for the anticipated thermal initial conditions with 9i — 0{1) p^ . 

As applied to the thermal scenario, the expression (||) is actually a 
very considerable underestimate for at least two reasons: First, the axions 
are not 'created' by the mass 'switch on' at t = tgco, they are already 
there with a specific momentum spectrum g{k) determined by dynamical 
mechanisms prior to this time. The actual axion number obtained from g{k) 
is much larger than the rms average assumed in (^) which ignores the true 
particle content. Secondly, this estimate was derived before much stronger 
topological effects were realized, notably the presence of axion strings and 
domain walls. In any case, these nonlinear effects complicate the oscillatory 
behaviour considerably, implying that the homogeneous estimate (^) is 
poorly motivated. 

2.2. Axion string network decay 

Axions and axion strings are inextricably intertwined. Like ordinary su- 
perconductors or superfluid ^He, axion models contain a broken U{1)- 
symmetry and so there exist vortex-line solutions. Combine this fact with 
the PQ phase transition, which means the field is uncorrelated beyond the 
horizon, and a random network of axion strings must inevitably form. An 
axion string corresponds to a non-trivial winding from to 2tt of the axion 
field 6 around the bottom of its 'Mexican hat' potential. It is a global 
string with long-range fields, so its energy per unit length /i has a logarith- 
mic divergence which is cut-off by the string curvature radius R'^t, that 
is, n « 27r/j^ ln(t/(5) , where the string core width is 5 w f~^. The axion 
string, despite this logarithmic divergence, is a strongly localized object 
and is likely to behave as such; if we have a string stretching across the 
horizon at the QCD temperature, then \n.{t/5) ~ 65 and over 95% of its 
energy lies within a tight cylinder enclosing only 0.1% of the horizon vol- 
ume. To first order, then, the string behaves like a local cosmic string, a 



cosmology) [|, |o|: 





3 



fact that can be established by a precise analytic derivation and careful 
comparison with numerical simulations . 

After formation and a short period of damped evolution, the axion 
string network will evolve towards a scale-invariant regime with a fixed 
number of strings crossing each horizon volume (for a cosmic string review 
see ref. This gradual demise of the network is achieved by the pro- 

duction of small loops which oscillate relativistically and radiate primarily 
into axions. The overall density of strings splits neatly into two distinct 
parts, long strings with length £ > t and a population of small loops i < t, 
that is, p = Poo + Pl- High resolution munerical simulations confirm this 
picture of string evolution and suggest that the long string density during 
the radiation era is Poc ~ 13/i/t^ To date, analytic descriptions of the 
loop distribution have used the well-known string 'one scale' model, which 
predicts a number density of loops defined as p£n{£,t)di ~ pL{£,t)M in 
the interval £ to £ -|- (i£ to be given by 

(£ + «t)5/2t3/2 ' 

where a is the typical loop creation size relative to the horizon and k w 
65/[27rln(i/(5)] is the loop radiation rate iQ. Once formed at i = with 
length £oi a typical loop shrinks linearly as it decays into axions i = — 
Kit — to). The key uncertainty in this treatment is the loop creation size a, 
but compelling heuristic arguments place it near the radiative backreaction 
scale, a ~ K. (If this is the case, we note that the loop contribution is over 
an order of magnitude larger than direct axion radiation from long strings.) 

String loops oscillate with a period T = £/2 and radiate into harmonics 
of this frequency (labelled by n), just like other classical sources. Unless 
a loop has a particularly degenerate trajectory, it will have a radiation 
spectrum P„ oc n^'^ with a spectral index q > 4/3, that is, the spectrum 
is dominated by the lowest available modes. Given the loop density (^), 
we can then calculate the spectral number density of axions dnn/duj, which 
turns out to be essentially independent of the exact loop radiation spectrum 
for q > 4/3. From this expression we can integrate over uj to find the total 
axion number at the time tqcD, that is, when the axion mass 'switches on' 
and the string network annihilates. Subsequently, the axion number will be 
conserved, so we can find the number-to-entropy ratio and project forward 
to the present day. Multiplying the present number density by the axion 
mass TTia yields the overall axion string contribution to the density of the 
universe @: 

/10~^eV\ ^'^^ 

f^a,stri„g « llOA/i-2 f{a/K) , (4) 
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where 

\ 3/2 



fia/K) 
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(5) 



The key additional uncertainty from the string model is the ratio aj k ^ 
0(1), which should be clearly distinguished from particle physics and cos- 
mological uncertainties inherent in A and h (which appear in all estimates 
of r^a)- With a Hubble parameter near h = 0.5, the string estimate (^ 
tends to favour a dark matter axion with a mass ma ~ 100/ieV, more- 
over a comparison with (^ confirms that f^a. string is well over an order of 
magnitude larger than the 'misalignment' contribution. 

2.3. Hybrid defect annihilation 

Near the QCD phase transition the axion acquires a mass and network evo- 
lution alters dramatically because domain walls form. Large field variations 
around the strings collapse into these domain walls, which subsequently be- 
gin to dominate over the string dynamics. This occurs when the wall sur- 
face tension a becomes comparable to the string tension due to the typical 
curvature a ^/t. The demise of the hybrid string-wall network proceeds 
rapidly, as demonstrated numerically The strings frequently intersect 
and intercommute with the walls, effectively 'slicing up' the network into 
small oscillating walls bounded by string loops. Multiple self-intersections 
will reduce these pieces in size until the strings dominate the dynamics 
again and decay continues through axion emission. 

An order-of-magnitude estimate of the demise of the string-domain wall 
network indicates that there is an additional contribution 

/in-6pV\ ^-^^ 
^^a,dw ~ O(10)A/i-2 p^^^^ . (6) 

This 'domain wall' contribution is ultimately due to loops which are cre- 
ated at the time ^ iqcD- Although the resulting loop density will be 
similar to (||), there is not the same accumulation from early times, so it 
is likely to be subdominant |14 relative to (^. More recent work, 



questions this picture by suggesting that the walls stretching between long 
strings dominate and will produce a contribution anywhere in the wide 
range Oa,dw ~ (l-44)r^a,string; however, this assertion requires stronger 
quantitative support. Overall, like most effects,^ the domain wall contri- 
bution will serve to further strengthen the string bound (|^) on the axion. 

^ We note briefly that it is also possible to weaken any axion mass bound through 
catastrophic entropy production between the QCD-scale and nucleosynthesis, that is, in 
the timescale range 10~^s^tent ^W~^s. Usually this involves the energy density of the 
universe becoming temporarily dominated by an exotic massive particle with a tuned 
decay timescale. 
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Up to this point we have only considered the simplest axion models 
with a unique vacuum iV = 1, so what happens when iV > 1? In this case, 
any strings present become attached to N domain walls at the QCD-scale. 
Such a network 'scales' rather than annihilates, and so it is cosmologically 
disastrous being incompatible (at the very least) with CMB isotropy. 



3. Theoretical challenges to nig, lOO^eV 

The conclusion that that Wa is most likely around 100//eV is an impor- 
tant one since the only experimental search for axions with a realistic level 
of sensitivity can only detect axions if Wa ~ 1 — 10/ieV. The possibility 
of detection by the current generation of experiments might still be pos- 
sible, either by pushing the theoretical and cosmological uncertainties of 
the calculation presented in the previous section to the limit, allowing for 
entropy production by a late decaying massive particle, or by allowing an 
inflationary epoch which prevents PQ symmetry restoration!^ However, 
the theoretical framework which we have presented in the previous section 
has been also challenged conceptually and quantitatively in the literature 
and here we shall comment critically on this work. 

The basis of our understanding global string evolution and radiation 
presented in the previous section is based on the assumption that global 
strings behave in an essentially similar way to local strings described by 
the an extension of the Nambu action, known as the Kalb-Ramond action. 
This kind of action is used frequently in superstring theory to describe the 
dynamics of strings coupled to a massless scalar field, and its efficacy for 
describing global string dynamics was established quantitatively in ref. . 
Since the global strings have long range forces and a logarithmically diver- 
gent energy density, one might wonder how this is possible. But as we 
have already pointed out over 95% of the energy lies within 0.1% of the 
volume for cosmological scales. For the point of view of this discussion we 
will separate the challenges in the literature into the conceptual jl^ and 
the quantitative fl^ . Both are just supported by numerical simulations, 
rather than the analytic arguments and numerical calculations used in our 
description. The main point being that use of numerical simulations of 
limited physical size and dynamic range can be hazardous when making 
extrapolations to cosmological scales where the strength of the relevant 
effects is very much diminished. 

^ We should note that as discussed in ref. 1^, even if the inflationary reheat temper- 
ature is below /a, it is still possible for the strings to be formed and for them to be the 
dominant source of axions. 
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3.1. Conceptual issues : the spectrum of radiation from strings 



While the issue of the spectrum of radiation from global strings is a rich and 
interesting theoretical subject, its effects on the cosmological axion density 
are much simpler than often realised. The main issue, as elucidated in 
ref. [p^ , relates to the power in high frequency modes. The frequency of 
the nth radiation mode for a loop is w„ — 47m/Z, where I is its length. If the 
fall-off of the power spectrum of radiation is a power law (P„ oc for large 
n, as suggested by analytic calculations using the Kalb-Ramond action ]20[| , 
then so long as 9 > 1 the calculation we presented in the previous section 
is still valid. This is because only in the divergent case q = 1 does the 
calculation of the axion density depend sensitively on the high frequency 



cut-off |17 . A wide range of analytic and numerical calculations support 
this issue, but since this issue has been discussed extensively in a large 
number of publications, we see no reason to discuss it further here. 



3.2. Quantitative issues : the scaling density of strings 

Recent simulations of a network of global strings on a 'large' (256'^) grid 
has suggested that the scaling density of strings ^ — poot^fJ. « 1, rather 
than ^ « 13 as used in our calculation. This clearly has a significant effect 
on the axion mass prediction, reducing it to around lOneV, since the axion 
density depends linearly on ^. The argument that ^ « 13 is based on the 
premise that the cosmological scale global strings behave in an essentially 
similar way to the Nambu strings for the reasons already discussed, and 
this may have to be modified to account for any alight deviations from 
this. However, in the simulations presented in ref. ||l^ there is no reason 
to believe that the strings will behave anything like Nambu strings because 
log{t/S) is at most 4, that is, in the same 0.1% of the volume of space, 
there is less than 50% of the energy. The long range forces between strings 
will be more than 10 times stronger than in the cosmological context and 
will substantially effect the dynamics and hence the scaling density. 

Furthermore, the dynamic range of the simulations is much smaller than 
it appears. The time of formation of the string network is t{ = 20ti, where 
ti is the start of the simulation, and the time of the end of the simulation 
is tc = 80ii. Hence, in real time t the dynamic range is tc/tf — 4, not 80 
as one might think, and the true dynamic range in conformal time is only 
Te/rf = (tc/tf)^^^ = 2. Therefore, the simulations are too small to simulate 
realistic cosmological behaviour and have very limited dynamic range. 
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